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65 4% 
钟 时 间 尺 度 , 利用 Vondrak_Cepek (V-C) 平 滑 对 两 
类 时 间 尺 度 有 效 组 合 , 改善 守 时 型 原子 钟 综合 时 间 
尺度 准确 度 和 稳定 度 指 标 . 首先 基于 最 小 误差 的 理 
论 基 础 ,对 氢 钟 频 漂 参 数 预测 进行 数学 建 模 ， 以 更 
准确 地 估计 该 参数 ; 然后 利用 二 次 指数 平滑 方法 产 
生 全 氧 钟 时 间 尺 度 ; 最 后 利用 V-C 平 滑 方法 将 两 类 
原子 钟 组 时 间 尺 度 进行 结合 , 改善 守 时 型 原子 钟 综 
合 时 间 尺 度 性 能 . 
ed md 
1. 以 最 小 误差 理论 为 指导 , 建立 氧 原子 钟 频率 
漂移 参数 预测 模型 , 改善 氧 原子 钟 钟 差 预 测 的 准确 
2. 提出 在 时 间 尺 度 层面 对 不 同类 型 原子 钟 数 
据 进 行 融合 的 概念 , 设计 V-C 滤 波 方法 , 优选 参数 ， 
并 利用 全 和 毛 钟 时 间 尺 度 的 一 阶 导数 改善 全 钨 钟 时 
间 尺 度 性 能 , 最 终 实 现 氧 钨 综合 时 间 尺 度 长 短期 稳 
Et 


生 能 的 提升 . 
2 理论 基础 

二 次 指数 平滑 钟 差 预测 方法 

传统 的 二 次 指数 平滑 方法 的 主要 步骤 是 钟 差 
预测 、 钟 速 预测 及 每 台 原子 钟 权重 的 确定 . LA 
时 间 t 的 钟 差 测量 数据 为 依据 ,对 下 一 时 刻 的 钟 差 
进行 估计 ， 
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7) 是 原子 钟 ; 在 (t+7) 时 刻 的 钟 差 估计 . 同时 , 根据 原 
子 钟 钟 差 的 数据 特点 , 对 应 于 n 台 原子 钟 可 得 到 mn 个 
主 钟 与 “ 纸 面 时 "的 差 . 最 后 , 对 上 述 所 得 的 mn 个 主 钟 
与 “ 纸 面 时 ”的 差 进行 合理 的 加 权 产 生 最 优 的 时 间 
尺度 ， 
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其 中 X;(t 十 7) 是 原子 钟 ;在 (t 十 7) 时 刻 的 钟 差 , wj 为 
第 i 台 原 子 钟 权重 ， Xij(t + 7T) 表 示 在 (t 十 7) 时 刻 原 
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子 钟 ? 和 和 7 的 钟 差 测量 值 . Gd 的 预测 ; 
原子 钟 钟 差 预测 准确 度 , 进 
能 . 
2.2 ”优化 的 二 次 指数 平滑 方法 模型 
利用 预测 误差 动态 估计 氢 原 子 钟 频 率 漂 移 参 
数 ， 相 比 传统 二 次 指数 平滑 钟 差 估 计 ， 增 加 了 预 
测 间 隔 参 数 . 在 该 时 间 段 内 ,根据 预测 误差 RMS 
(Root Mean Square) 最 小 原则 对 4d 和 平均 时 间 进 行 
最 优 估计 , 从 而 提升 钟 差 估 计 的 准确 度 . 
预测 算 子 表达 式 如 下 


= a(t) + Tiy (t, T2) + iar? (1 + =) , 
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子 钟 钟 速 因 末 移 动 平 均 算 子 
预测 误差 定义 为 
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间隔 为 ni 时 的 钟 差 值 和 预测 值 
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影响 综合 时 间 尺 度 性 
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其中 入 (i) 表示 测量 
之 差 . 
预测 误差 RMS 为 
vs{eton m= | See(f) x|H(J2mf)|"df , 
0 
(5) 
其 中 了 为 原子 钟 频率 ,s,s 表示 原子 钟 钟 差 的 离 差 平 
方 和 . 传递 函数 | 瑟 (727 力 | 的 定义 如 下 ， 
IHGA) =4 (2) (1+ 2) sin® (rfr) + 
4 (1 + =) sin? (m fri) — 
4 (=) sin? (mf (r+71)) . 
2.3 Vondrak_Cepek 平 滑 建 模 
优化 的 三 次 指数 平滑 方法 产生 氧 钟 时 间 尺 度 ， 
改善 的 加 权 平 均 方法 得 到 全 钨 钟 时 间 尺 度 , 最 后 对 
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两 种 不 同类 型 的 时 间 尺 度 进行 有 机 结合 . 这 里 全 钨 
钟 时 间 尺 度 (TA-Cs) 用 m(t) 表 示 , 它 是 时 间 t 的 函数 ， 
在 t 时 刻 的 一 阶 导 数 定义 为 


m'(t) = [m (t+ Tes) — m(t)] /Tes ; 


种 新 的 守 时 型 原子 钟 综合 时 间 尺 度 方法 有 


多 1 2 


后 将 扣 去 钟 速 的 钨 钟 和 扣 去 钟 速 频 漂 的 氧 钟 组 成 
守 时 钟 组 , 利用 加 权 平 均 方法 产生 综合 时 间 尺 度 . 
] 氧 钟 作为 测量 基准 , 利 
计 氢 原子 钟 参数 产生 氧 


而 本 文 所 提 方 法 避免 了 使 
钟 组 时 间 尺 度 , 并 利用 加 权 平 均 算 法 产生 钨 钟 时 间 


其 中 7cs 表 示 钨 原子 钟 的 钟 差 测 量 间 隔 . MORR 
全 氧 钟 时 间 尺 度 (TA-H), 其 一 阶 导数 定义 为 


M'(t) = [M (t + ma) — M(H] /ma, 


其 中 表示 氧 原子 钟 的 钟 差 测量 间隔 . m (EFM (t) 
表示 两 类 钟 组 时 间 尺 度 速率 . 如 果 饮 原 子 钟 和 氧 原 
子 钟 测量 间隔 相同 , Wm (t) = M(t). 则 有 


MJD2 MJD2 
se 


M'((TA-H))d¢, 

(7) 
其 中 MJD (Modified Julian Day) 为 简化 儒 略 日 . 钨 
钟 组 时 间 尺 度 的 估计 计算 如 下 : 


m(MJD2) = m(MJD1) + M'(t)At, 


m'(TA-Cs)dt = f 


MJD1 


(8) 


其 中 At 表示 时 间 间 隔 . 

本 文 在 计算 时 , 选择 的 原子 钟 数据 测量 间隔 为 
1h, 因此 两 种 原子 钟 的 测量 间隔 相同 . AS Ze 
时 间 尺 度 在 下 一 时 刻 的 钟 差 估 计 可 以 使 用 此 时 全 
钨 钟 时 间 斥 度 钟 差 值 与 全 氧 钟 时 间 尺 度 钟 速 和 时 
间 间 隔 乘积 的 和 产生 . 原子 钟 数据 测量 间隔 的 选取 
也 对 时 间 尺 度 的 性 能 有 着 很 大 的 影响 . 测量 间隔 过 
K, 对 于 和 氢 原子 钟 来 说 在 该 测量 区 间 的 频 漂 变 大 ， 


用 改进 的 二 次 平滑 动态 估 

尺度 , 最 后 采用 V-C 平 滑 结合 两 类 钟 组 时 间 尺 度 的 
长 短期 优势 . 该 方法 利用 最 优 估计 方法 动态 预测 氢 
钟 频 漂 参 数 , 较 二 次 拟 合 方法 的 频 漂 参数 估计 更 为 
准确 合理 . 同时 , 利用 加 权 平 均 方 法 产生 钨 钟 时 间 
尺度 , 利用 二 次 平滑 方法 产生 氧 钟 时 间 尺 度 . 算法 
选择 时 考虑 了 加 权 平 均 方法 依据 原子 钟 所 有 历史 
信息 来 分 配 权重 , 更 注重 时 间 尺 度 的 长 期 性 能 , 适 
合 长 期 性 能 更 好 的 钨 原子 钟 组 ; 二 次 指数 平滑 方法 
则 更 关注 时 间 尺 度 实 时 性 , 更 适合 短期 性 能 良好 的 
ARTH. 最 终 采 用 V-C 联 合 平滑 方法 综合 两 类 
时 间 尺 度 优势 , 提升 综合 时 间 尺 度 性 能 . 


3 ”实验 及 结果 分 析 
3.1 ”算法 设计 

本 文 基于 Fortran 语 言 , 在 Visual Studio 编 译 环 
境 开发 了 一 套 软件 系统 来 实现 改进 的 二 次 指数 平 
滑 (OX) 和 V-C 平 滑 (OX+V-C) 氢 钨 综合 时 间 尺 度 
算法 . 图 1 为 OX 十 V-C 方 法 设计 过 程 ，OX+V-C 实 
现 主要 包括 如 下 几 个 环节 : q 的 估计 、n 台 原子 钟 
与 “ 纸 面 时 ”的 时 差 预 测 、 每 台 原 子 钟 取 权 、 主 钟 
与 “ 纸 面 时 ”的 时 差 预测 及 钟 速 预测 、 氧 钟 组 改进 
二 次 指数 平滑 时 间 尺 度 、 钨 钟 组 加 权 平 均 时 间 尺 
度 、 全 氧 钟 时 间 尺 度 钟 速 产 生 、V-C 滤 波 器 设计 及 


综合 时 间 尺 度 准 确 性 将 受 影响 ; 测量 频率 过 大 , 原 
子 钟 自身 的 噪声 可 能 被 测量 噪声 所 掩盖 , 导致 原子 
钟 准确 建 模 更 加 困难 , 最 终 影响 时 间 尺 度 性 能 . 
2.4 ”算法 优势 

文献 8] 中 产生 氧 钨 综合 时 间 尺 度 的 方法 是 以 
高 性 能 氧 原子 钟 作 为 测量 参考 , 获取 每 台 钨 原子 钟 
相对 于 该 高 性 能 氨 原 子 钟 的 钟 差 , 然后 对 每 台 钨 原 
子 钟 钟 差 数 据 利 用 数学 滤波 方法 去 除 噪声 , 再 利用 
最 小 二 乘法 对 钨 钟 钟 差 数 据 进 行 拟 合 , 估计 每 台 钨 
原子 钟 钟 速 , 并 将 其 扣除 以 备 使 用 . TREY SU 


t 


结果 评估 . 
3.2 ”实验 过 程 及 结果 分 析 
本 实验 以 氧 原子 钟 和 钨 原子 钟 为 研究 对 象 ， 
3 台 氢 钟 和 3 台 钨 钟 构成 守 时 钟 组 . 选择 MJD59945- 
59975 (2023 年 1 月 1 日 到 1 月 31 日 ) 的 原子 钟 数据 进 
行 本 文 所 提 方 法 有 效 性 验证 . 首先 , 根据 最 小 误差 
理论 估计 每 台 氧 原子 钟 频率 漂移 参数 并 设计 指数 
滤波 器 ,计算 改进 二 次 指数 平滑 时 间 尺 度 TA-DH. 
为 分 析 该 方法 对 时 间 尺 度 性 能 的 提升 , 与 传统 二 
次 指数 平滑 的 全 氧 钟 时 间 尺 度 进行 了 对 比 . 根据 


频率 和 频 漂 参 数 是 利用 二 次 拟 合 方法 进行 预测 , 最 
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氨 原 子 钟 最 后 5 d 钟 差 数 据 的 二 次 拟 合 对 频率 漂移 
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参数 进行 估计 , 并 计算 二 次 指数 平滑 时 间 尺 度 TA- (ADEV), 结果 如 表 1 所 示 . 
H, 结果 如 图 2. 同时 计算 两 种 时 间 尺 度 的 阿 伦 偏差 
Cs clock1 
IWAM 
Cs clock2 > TA-Cs 
Cs clock3 TA- Result 
> SFO HCs 7] evaluation 
S 
Hm 1 
EEN Frequency 
Hm 2 TA-DH ~> of TA- 
DH 
Hm 3 
41 OX+V-CA HR AI ERRE RATE 
Fig.1 The algorithm flow of the comprehensive time scale of OX+V-C hydrogen masers and cesium atomic clocks 
Note: Hm in the figure represents the hydrogen master, TA-DH represents the time scale of the hydrogen clock generated 


after dynamically estimating the frequency drift parameters, IWAM represents the improved weighted average method, IEFM 


represents the improved quadratic smoothing method, SFO represents the optimal smoothing factor, TA represents timescale, 


and TA-HCs represents the comprehensive time scale of hydrogen maser and cesium atomic clock. 


TA-DH 
TA-H 


Time difference/ns 


-8 f f f f f 1 
59945 59950 59955 59960 59965 59970 59975 59980 
MJD 


次 指数 平滑 时 间 尺 度 和 二 次 指数 平滑 时 间 


2 


43 


R 


Fig.2 Improved quadratic exponential smoothing time scale 


and quadratic exponential smoothing time scale 


由 图 2 可 知 , 改进 的 二 次 指数 平滑 得 到 的 时 间 
尺度 在 31 d 的 时 间 间 隔 中 ， 偶 差 在 4 ns 以 内 ， 而 固 
定 频 漂 方 法 产生 的 时 间 尺 度 在 相同 时 间 段 内 的 偏 
差 在 8 ns 之 内 . 故 所 提 改 进 的 二 次 指数 平滑 方法 
可 有 效 改善 全 氧 钟 时 间 尺 度 的 准确 度 . 由 表 1 得 ， 
TA-DH 较 TA-H 在 不 同 的 平均 时 间 下 的 阿 伦 偏差 值 


E. TA-DH 在 1 h 的 稳定 度 为 1.29 x 1071, 天 稳定 
度 为 7.0 x 10715. 综 上 所 述 , 改进 的 二 次 指数 平滑 方 
法 可 以 有 效 提 升 时 间 尺 度 的 准确 度 , 同时 对 稳定 度 
指标 也 有 所 改善 . 

表 1 改进 二 次 指数 平滑 时 间 尺 度 和 二 次 指数 平滑 时 间 尺 度 


阿 伦 偏差 
Table 1 The Allan deviation (ADEV) of 


improved quadratic exponential smoothing time 


scale and quadratic exponential smoothing time 


scale 
ADEV 
Averaging time/ns 
TA-H TA-DH 
3.60 x 1077 1.30 x 10714 1.29 x 10714 
7.20 x 1012 7.54x 107 667x107! 
1.44 x 1018 3.50 x 10715 3.38 x 107% 
2.88 x 101% 3.82 x 107! 3.75 x 1071° 
5.76 x 1013 6.90 x 1071+ 6.90 x 1071° 
1.15 x 1014 9.96 x 10715 9.67 x 107% 
2.30 x 1014 9.60 x 107! 9.46 x 10715 
4.61 x 1014 7.08 x 10715 6.87 x 107% 


小 ， 说 明 改 进 的 二 次 指数 平滑 时 间 尺 度 长 短期 稳 
定 度 略 优 于 相同 钟 组 采用 传统 方法 计算 的 时 间 尺 
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然后 , 通过 加 权 平 均 算法 产生 钨 钟 组 时 间 尺 度 


BS 


= HWH: 


TA-Cs. 全 钨 钟 时 间 尺 度 短 期 波动 大 ,而 长 期 波动 
较 小 ; 全 氧 钟 时 间 尺 度 则 相反 , 这 种 情况 与 V-C 平 
滑 理论 模型 非常 相符 . 最 后 , 通过 计算 全 氧 钟 时 间 
尺度 一 阶 差分 , 根据 V-C 平 滑 算 法 计算 氧 饮 综 合 时 
间 尺 度 , 时 间 尺 度 曲线 见 图 3. 图 4 是 3 台 钨 钟 加 权 平 
均 产生 的 时 间 尺 度 与 OX+V-C 时 间 尺 度 残 差分 布 
情况 .从 图 3 可 得 , 两 条 曲线 走向 一 致 ， 但 OX 十 V- 
C 曲 线 较 TA-Cs 平 滑 很 多 , 证 明 该 方法 有 效 改善 了 
时 间 尺 度 的 短期 波动 且 保 留 了 全 饮 钟 时 间 尺 度 的 
长 期 性 能 . 同时 ， 采用 全 钨 钟 时 间 尺 度 与 OX 十 V- 
C 时 间 尺 度 之 间 的 均 方 根 误 差 (Root Mean Square 
Error，RMSE) 对 实验 结果 进行 评估 ，RMSE = 


0.02， 表 明 V-C 平 滑 方法 得 到 的 时 间 尺 度 能 够 跟 


种 新 的 守 时 型 原子 钟 综合 时 间 斥 度 方 法 和 


vi 1 


,曲线 代表 AT1 方 法 产生 的 氧 钨 综合 时 间 尺 度 , 绿 
和 线 表示 本 文 所 提 方 法 OX+TV-C 产 生 的 氧 钨 综 
间 尺 度 ， 黑 色 曲 线 为 Kalman 滤 波 算法 产生 的 
AiG SRA IN TARE. 从 图 5 可 看 出 ，OX+V-C 时 间 
尺度 曲线 比 其 他 3 条 平滑 , 且 在 一 个 月 的 时 间 间 隔 
内 ， 时 间 偏 差 量 最 小 . 同时 ， 计 算 以 上 4 种 综合 时 
间 尺 度 的 阿 伦 偏差 , 数据 指标 和 分 布 情况 见 表 2 和 
图 6. 从 表 2 和 图 6 可 得 ，OX 十 V-C 时 间 尺 度 在 1 h 稳 
定 度 为 1.60 x 10-5， 天 稳定 度 为 3.0 x 1071, 
于 AT1、ALGOS、Kalman 滤 波 3 种 方法 产生 的 时 
间 尺 度 对 应 指标 . 因此 , 进一步 说 明 论 文中 所 提 方 
法 的 长 短期 稳定 度 均 优 于 3 种 经 典 方 法 得 到 的 时 间 


随 TA-Cs 的 趋势 , 保留 其 大 部 分 信息 ， 以 有 效 利用 
全 钨 钟 时 间 尺 度 长 期 优势 . 图 4 可 得 , 两 种 时 间 尺 度 
残 差 在 0 值 左右 波动 ,其 残 差 也 服从 正 态 时 间 分 布 
规律 , 这 与 白 噪 声 特 性 吻合 . 


La 
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Fig.3 The comprehensive time scale of hydrogen masers and 


cesium atomic clocks and time scale of all cesium clocks 


为 检验 OX+V-C 氢 钨 综合 时 间 尺 度 方法 长 短 
期 稳定 度 是 否 有 所 改善 , 本 文 也 基于 AT1、ALGOS 
以 及 Kalman 滤 波 的 3 种 经 典 算法 产生 的 综合 时 间 
尺度 相应 指标 进行 了 比较 . 首先 将 3 台 氧 钟 扣 去 频 
率 漂 移 对 钟 差 的 影响 ,再 与 3 台 钨 钟 参与 计算 ， 并 
利用 传统 ATI1、ALGOS 和 Kalman 滤 波 3 种 算法 产 
生 的 综合 时 间 尺 度 . 时 间 尺 度 结果 见 图 5. 红色 曲 
线 代 表 ALGOS 方 法 产生 的 氧 饮 综 合 时 间 尺 度 ， 蓝 


尺度 稳定 度 . 同时 证 明了 该 方法 能 有 效 利 用 两 类 原 
子 钟 组 时 间 尺 度 优势 , 进一步 提升 综合 时 间 尺 度 性 
能 . 


0 
-0.5 0 0.5 
Residual/ns 


尺度 和 全 饮 钟 时 间 


图 4， SERATA 


尺度 残 差 分 布 


Fig.4 Time scale residual distribution of hydrogen maser 


and cesium atomic clock and all cesium clock 


同时 , 本 文 也 利用 同样 的 守 时 钟 组 原子 钟 数据 ， 
依据 文献 81 和 文献 9 中 采用 的 氧 钨 综合 时 间 尺 度 
计算 方法 分 别 产生 两 种 融合 时 间 尺 度 TA1 和 TA2; 
利用 优化 的 二 次 平滑 方法 分 别 产生 氧 钟 组 和 钨 钟 
组 时 间 尺 度 ， 再 经 V-C 平 滑 产 生 综 合 时 间 尺 度 TA- 
EF-VC. 分 别 计算 3 类 时 间 尺 度 的 稳定 性 指标 ， 并 
与 本 文 所 提 方 法 进行 了 比较 , 结果 见 表 3. 从 表 3 中 
可 看 出 , 本 文 所 提 方 法 OX+V-C 产 生 的 综合 时 间 尺 
度 时 稳 及 天 稳 都 高 于 其 他 3 种 时 间 尺 度 的 对 应 指标 ， 
进一步 证 明了 OX+V-C 算 法 在 本 文选 定 的 6 台 原 子 
钟 数据 计算 中 性 能 最 优 . 而 本 文 所 提 方 法 产生 的 时 
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间 尺 度 稳定 性 指标 也 比 TA-EF-VC 方 法 更 好 , 分 析 法 ,建立 一 种 新 的 守 时 型 原子 钟 时 间 尺 度 产生 模 
原因 是 由 于 钨 钟 无 频 漂 项 , 优化 的 二 次 平滑 方法 主 型 , 通过 氧 钟 组 时 间 尺 度 的 导数 信息 改善 饮 钟 组 
要 针对 氢 钟 的 频 漂 进行 最 优 估 计 , 关注 时 间 尺 度 的 时 间 尺 度 短期 噪声 大 的 问题 ， 实现 了 全 氢 钟 和 全 


H 


短期 稳定 度 , 更 适合 短期 性 能 较 好 的 氧 钟 组 . 钨 钟 在 时 间 尺 度 层面 的 结合 . 实验 结果 表明 , 本 文 
所 提 方 法 产生 的 综合 时 间 尺 度 1 h 稳 定 度 达 1.60 x 
10- 革 ,天 稳定 度 达 3.0 x 10715, 较 ALGOS、AT1 


4+ we Se ` 
和 Kalman 滤 波 3 种 传统 方法 计算 的 综合 时 间 尺 度 
a 3f 长 短期 性 能 有 明显 提高 , 达到 了 预期 效果 , 证 明了 
3 2} 本 文 所 提 方 法 的 有 效 性 . 
2 
= ? 
E 0 10 7 
请 1l 
一 一 ALGOS 
| 一 一 AT 
-2 —— OX+V-C 
3 f i l l ， Kalman me 
59945 59950 59955 59960 59965 59970 59975 59980 内 104 F 
MJD < 
图 5 ALGOS、AT1、Kalman 和 OX 二 V-C 算 法 的 氨 饮 综合 时 间 尺 
度 
Fig.5 Comprehensive time scales of hydrogen maser and 10 o 1013 ol 1015 


cesium atomic clock for ALGOS, AT1, Kalman, and As . iie 
OX+V-C methods VEraging timers 
图 6 ALGOS、AT1、Kalman 和 OX 十 V-C 算 法 方法 综合 时 间 尺 度 
` Bil $63 Ah Ae R 
4 结论 
Fig.6 The trend of Allan deviation of comprehensive time 


本 文 基于 优化 的 二 次 指数 平滑 和 V-C 平 滑 方 scale of ALGOS, AT1, Kalman, and OX+V-C methods 


#2 ALGOS、AT1、Kalman 和 OX 十 V-C 算 法 方法 综合 时 间 尺 度 稳定 度 指标 
Table 2 The stability indicator of comprehensive time scale for ALGOS, AT1, Kalman, and OX+V-C 


methods 
ADEV 
Averaging time/ns 
ALGOS AT1 Kalman OX+V-C 
3.60 x 1012 5.22 x 10714 1.10 x 10714 2.31x107'* 1.60 x 1071 
7.20 x 101? 3.76 x 10714 5.58 x 10715 1.73x10- 1.39 x 10715 
1.44 x 101 2.77 x 107! 3.29 x 10715 8.65x 10715 2.14 x 10715 
2.88 x 101% 1.83 x 1074 3.35 x 10715 6.34 x 10715 2.39 x 1071* 
5.76 x 101% 1.11 x 10714 408x107 7.61x107'° 3.03 x 107% 
1.15 x 1014 4.03 x 10715 4.49 x 10715 4.33 x 10715 3.46 x 1071 
2.30 x 1014 3.32x 107 3.55x107% 347x10 ™ 3.09x 107% 
4.61 x 104 3.73 x 107! 4.06 x 10715 4.03 x 10715 2.40 x 10715 
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#3 TA1, TA2, TA-EF-VC, OX+V-C 
Table 3 The stability indicator of comprehensive time scale of TA1, TA2, TA-EF-VC and OX+V-C 


萌 : 


种 新 的 守 时 型 原子 钟 综合 时 间 尺 度 方法 有 


= 
综合 


究 


时 间 尺 度 稳定 度 指标 


Averaging time/ns SaN 
TA1 TA2 TA-EF-VC OX+V-C 
3.60 x 10? 8.43 x 10714 4.32 x 10714 3.31 x 10714 1.60 x 1071 
7.20 x 10"? 6.38 x 107'* 3.58 x 107'4 2.33 x 107" 1.39 x 107" 
1.44 x 1013 2.93 x 10714 2.23x107'4 6.72x 10715 2.14 x 107" 
2.88 x 10" 2.36 x 10714 1.95 x 10714 5.23 x 10715 2.39 x 107" 
5.76 x 1013 1.92 x 1071 1.06 x 10714 8.84x 1071" 3.03 x 107" 
1.15 x 1014 1.32 x 10714 9.89 x 10715 5.62 x 10715 3.46 x 1071 
2.30 x 10" 9.37 x 10715 7.35 x 1071 4.48 x 10715 3.09 x 107" 
4.61 x 1014 6.64 x 10715 5.43 x 1071" 3.29 x 10715 2.40 x 107" 
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A New Method of Comprehensive Time Scale for Time-keeping 
Atomic Clock 
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ABSTRACT A comprehensive time scale calculation method for hydrogen maser and cesium atomic clock 
based on the improved exponential filtering and Vondrak_Cepek joint smoothing method is proposed. It is 
based on the minimum error method, which dynamically estimates the frequency drift parameters of the 
hydrogen maser and improves the accuracy of predicting the clock bias of the hydrogen atomic clock; With 
the method of improved quadratic exponential smoothing, the hydrogen masers time scale is generated, 
and the weighted average method is used to generate the cesium atomic clocks’ time scale. At the same 
time, the Vondrak_Cepek filter is designed to combine the long-term and short-term stability advantages 
of the two types of time scales to improve the comprehensive time scale performance. Experimental results 
show that the proposed method produces hydrogen cesium with a comprehensive time scale stability of 
up to 1.60 x 10715 h~!, stable up to 3 x10715 d~', outperforming the time-scale performance generated 
by three classic methods: ALGOS, AT1, and Kalman filtering. 


Key words time scale, clock difference prediction, dynamic estimation, stability 
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